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Abstract 

We present the algebraie representation and basie algorithms for Multi- 
Aspeet Graphs (MAGs). A MAG is a strueture eapable of representing mul¬ 
tilayer and time-varying networks, as well as higher-order networks, while 
also having the property of being isomorphie to a direeted graph. In par- 
tieular, we show that, as a eonsequenee of the properties assoeiated with 
the MAG strueture, a MAG ean be represented in matrix form. Moreover, 
we also show that any possible MAG funetion (algorithm) ean be obtained 
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from this matrix-based representation. This is an important theoretical result 
since it paves the way for adapting well-known graph algorithms for appli¬ 
cation in MAGs. We present a set of basic MAG algorithms, constructed 
from well-known graph algorithms, such as degree computing. Breadth First 
Search (BFS), and Depth First Search (DFS). These algorithms adapted to 
the MAG context can be used as primitives for building other more sophisti¬ 
cated MAG algorithms. Therefore, such examples can be seen as guidelines 
on how to properly derive MAG algorithms from basic algorithms on di¬ 
rected graph. We also make available Python implementations of all the 
algorithms presented in this paper. 


1 Introduction 


Graph theory finds many applications in the representation and analysis of com¬ 
plex networked systems [[T]-j^. In most cases, the utility of the graph abstraction 
comes from its inherent ability to represent binary transitive relations (i.e. tran¬ 
sitive relations between two objects), which due to the transitivity property gives 
raise to key concepts, such as walks, paths, and connectivity. This graph concep¬ 
tual framework allowed the emergence of basic algorithms, such as Breadth First 
Search (BFS) and Depth First Search (DFS) [Ql^. These basic graph algorithms, 
in their turn, made possible the development of more sophisticated algorithms 
for the analysis of specific properties of complex networks, such as network cen¬ 
trality or network robustness [[6]-Q, and also the analysis of dynamic processes 
in complex networks, such as network generative processes or information diffu¬ 
sion (MU- Several generalizations of the basic graph concept have been pro¬ 
posed for modelling complex systems that can be represented by layers of distinct 
networks p?|[T^ and also complex systems in which the network itself evolves 
with time [17,18 j. 


In our previous work [T^, we formalize the MultiAspect Graph (MAG) struc¬ 
ture, while also stating and proving its main properties. The adopted adjacency 
concept in MAGs is similar to the one found in simple directed graphs, where 
the adjacency is expressed between two vertices, leading to a structure in which 
an edge represents a binary relation between two composite objects. Moreover, 
in 10. we show that MAGs are closely related to simple directed graphs, as we 
prove that each MAG has a simple directed graph, which is isomorphic to it. This 
isomorphism relation between MAGs and directed graphs is a consequence of the 
fact that both MAGs and directed graphs share a similar adjacency relation. 
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MAGs find application in the representation and analysis of dynamic complex 
networks, sueh as multilayer or time-varying networks; or even networks that are 
both multilayer and time-varying as well as higher-order networks [ 20pT |. Exam¬ 
ples of sueh networks inelude faee-to-faee in-person eontaet networks [22|, mo¬ 
bile phone networks [[2^[2^, gene regulatory networks [|25||, urban transportation 


networks p6| , brain networks p7|[28| , soeial networks p9| , among many oth¬ 
ers. In particular, we have previously applied the MAG abstraetion from [1^ to 
different purposes, sueh as modeling time-varying graphs [j^, studying time een- 
trality in dynamic complex networks [^, and investigating social events based 
on mobile phone networks To illustrate the MAG eoneept in more details in 
this paper, we present in Section an example of modeling a simple illustrative 
multimodal urban transportation network. 

In this paper, we build upon the basie MAG properties presented in [19| and 
show that MAGs ean be represented by matriees in a form similar to those used 
for simple direeted graphs (i.e., those with no multiple edges). Moreover, we here 
show that any algorithm (funetion) on a MAG ean be obtained from its matrix 
representation. This is an important theoretieal result sinee it paves the way for 
adapting well-known graph algorithms for applieation in MAGs, thus easing the 
effort to develop the analysis and application of MAGs for modelling eomplex 
networked systems. We then present the most eornmon matrix representations 
that ean be applied to MAGs, although we do not detail all the properties of these 
matriees, sinee they are well established in the literature [|5 33|. Further, we in- 
troduee in detail the eonstruetion of MAG algorithms for eomputing degree, BPS, 
and DFS to exemplify how MAG algorithms ean be derived from traditional graph 
algorithms, thus providing an illustrative guideline for developing other more so- 
phistieated MAG algorithms in a similar way. As a further eontribution, we also 
make available Python implementations of all the algorithms presented in this 
paper at the following URF: http : / /github . com/wehmuthklaus/MAG_ 
Algorithms, 

This paper is organized as follows. Section [^briefly presents the basic MAG 
definitions and properties derived from p9| in order to allow enough baekground 
of the current paper. Seetion also presents illustrative examples of MAGs and 
its adjaeeney notion. Seetionshows the representation of MAGs by means of al- 
gebraie struetures, sueh as matriees. Emphasis is given to matrix representations, 
whieh are derived from the isomorphism relation between MAGs and simple di¬ 
reeted graphs. In partieular, we also introduee in Seetion [3T] the eompanion tuple, 
whieh is a eomplement to the MAG matrix representations. In Seetion Q we 
present basie MAG algorithms which are derived from well-known simple graph 
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algorithms. Further, in Section 4.2 we show that any algorithm (function) that 
can be defined for a MAG can be also obtained from its adjacency matrix and 
companion tuple, establishing the theoretical basis for deriving MAG algorithms 
from well-known simple graph algorithms. Finally, Section [^presents our final 
remarks and perspectives for future work. 


2 MultiAspect Graph (MAG) 


In this section, we present a formal definition of a MAG, as well as some key 
properties, which are formally stated and proved in p9|. 


2.1 MAG definition 

We define a MAG as H = (A^E), where £ is a set of edges and A is a finite list 
of aspects. Each aspect (p s A is a finite set, and the number of aspects p = \A\ is 
called the order of H. Each edge e e E is a tuple with 2 x p elements. All edges 
are constructed so that they are of the form (ai,... ,,..., bp), where ai , bi 
are elements of the first aspect of H, 02,^2 are elements of the second aspect 
of H, and so on, until ap,bp which are elements of the p-th aspect of H. Note 
that the ordered tuple that represents each MAG edge is constructed so that their 
elements are divided into two distinct groups, each having exactly one element of 
each aspect, in the same order as the aspects are defined on the list A. 

As a matter of notation, we say that A{H) is the aspect list of H and E{H) is 
the edge set of El. Eurther, A{H)[n] is the n-th aspect in A{H), \A{H)\n\ | = T„ is 
the number of elements in A(//)[n], and p = \A{H) \ is the order of H. 

In addition to the former definition, we define the following two sets con¬ 
structed from the cartesian products of aspects of an order p MAG: 

p 

¥(//)= ( 1 ) 

n=l 

the cartesian product of all the aspects of the MAG H, and 

2p 

E{H) = A{H)[{n —l){mod p) + i], (2) 

which is the set of all possible edges in the MAG H, so that E{H) c E(//). 
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We call u e Y{H) a composite vertex of MAG H. As a matter of notation, 
a composite vertex is always represented as a bold lowercase letter, as in u, for 
instance. From the properties stated for the MAG edge in our definition, it follows 
that an MAG edge is closely related to an ordered pair of composite vertices. For 
any given MAG H, every MAG edge esE (H) has the form (ai,..., ,..., bp), 

so that (ai,..., Op) e V (//) and bp) e V [H) are composite vertices of this 

given MAG H. From this, we can define two functions 

Tlo : E(//) ^ Y{H) (3) 

e = {ai,a2,...,ap,bi,b2,...,bp) ^ {01,02,... ,ap) = u, 

and 

(4) 

e = {oi,02,...,Op,bi,b2,...,bp) ^ {bi,b2,... ,bp) = v. 

We call Ttoie) the origin composite vertex of e and 7rj(e) the destination composite 
vertex of e. Moreover, we can define the function 

V/:E(//)-V(//) X V(//) (5) 

{no{e),nd{e)) = {{oi,... ,Op),{bi,... ,bp)) = (u,v), 

from which we can construct a directed graph Gh = {V{H),\i/{E{H)). In , we 
demonstrate that the directed graph Gh = {V{H), ^f{E{H)) is isomorphic to the 
MAG El from which it was originated. At this point, we can therefore define the 
function 


g : {A{H),E{H)) - {\{H),Y{H) X V(//)) (6) 

H^{Y{H),^f{E{H)), 

which maps every MAG H to its isomorphic directed graph g{H). Further, we 
define the set of functions 


%r.^{H)^A{H)[i\ 

{o\,02, ■ ■ ■, Op) I >■ Oi, 

which extracts the n-th element of a composite vertex tuple. 


( 7 ) 
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2.2 MAG sub-determination 

The sub-determination is a generalization of the aggregation concept applied to 
multilayer or time-varying graphs, in which all layers can be aggregated, result¬ 
ing in a traditional graph. Since a MAG can have more than 2 aspects, the sub¬ 
determination can be done in more ways than the aggregation. 

A given MAG H of order p, can be sub-determined in 2^ — 2 ways. For each of 
these IP — 2 ways, we have a list Ac(//) of the aspects used to determine 

an equivalence class. Note that in a MAG of order p = I (i.e. a traditional graph), 
a vertex can not be sub-determined, since 2P — 2 = 0. 

2.2.1 Sub-determined composite vertices 

Let with 1 < ^ < 2^ — 2, be an index for one of the possible ways to construct 
a proper nonempty sublist of aspects. From this, we can define a canonical rep¬ 
resentation of the sub-determination directly defined by For any given we 
consider the p-bit binary expansion of ^ that is used as an indicator showing which 
aspects of the original MAG are present on the sub-determination. More specifi¬ 
cally, the least significant bit indicates the presence or absence of the first aspect 
and the most significant bit indicates the presence or absence of the last aspect. 
By this convention, in a MAG with p = 3 aspects, we have that ^ = 0 0 12 corre¬ 
sponds to the sub-determination where only the first aspect is present, ^ = 01 02 
corresponds to the sub-determination where only the second aspect is present, 
^ = 10 12 corresponds to the sub-determination where both the first and the third 
aspects are present, and so on. By using this convention, we can directly associate 
a given ^ to its corresponding aspect sublist. 

Therefore, for each we have a unique sublist of aspects, such that 

Pi^ = \Ai^{H) \ is the order of the sub-determination We now define the set 

V((H)- X >!{(«)["]. (8) 

n=\ 

where Ni^{H) is the cartesian product of all the aspects in the sublist Ai^{H) of 
aspects, according to the index We call e a sub-determined vertex, 

according to the sub-determination 
We can now define the function 

5^ : V(//) ^ V^(//) (9) 

(ai,a2,...,ap) 
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where m = pi^. maps a composite vertex u e \{H) to the corresponding sub¬ 
determined composite vertex e \according to the sub-determination 
As eV^(//), it follows that a eA^(//)[l],... eA^(//)[m]. 

From the definition, it can be seen that the function 5,^ is not injective. Hence, the 
function for a given sub-determination can be used to define a equivalence 
relation in V(//), where for any given composite vertices u,v e V(//), we 
have that u v if and only if 5^(u) = S^(v). 

2.2.2 Sub-determined edges 

From the sub-determination ^ of order we can also construct the set 

2xp^ 

= 'X Ai^{H)[{n-l){mod pi^) + l], (10) 

n=l 

where = |A^(//)| is the order of the sub-determination and is the set 

of all possible sub-determined edges according to We then define the function 

E^:E{H) (11) 

where m = p^ and e (//) [1], ^P],..., e {H) [m] 

This function takes an edge to its sub-determined form according to ^ in a similar 
way as defined above for composite vertices. In general, the function is not in¬ 
jective. Consider two distinct edges ei,e 2 ^ E{H), such that ei and 62 differ only 
in aspects which are not in Ai^{H). Since Ei^{-) only contains values for aspects 
present in Ai^{H), it follows that E^{ei) = £^(^ 2 ), and therefore E^^ is not injec¬ 
tive. Further, consider an edge e e E (H) and its sub-determined edge = Ei^{e), 
such that TToiei^) = i-e. is a self-loop. Since self-loops are not allowed 

to be present on a MAG, it follows that ^ Ei^{E{H)). As consequence, we have 
thnt\E^iEiH))\^\EiH)\. 

2.2.3 Sub-determined MAGs 

For a given sub-determination ^ we have the sublist A ^(//) of considered aspects 
and also the sub-determined edges obtained from Based on them, we can now 
obtain a sub-determined MAG. For a given sub-determination ^ we define the 
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function 


: (A(//),£(//)) - (A^(//),E^(//)) (12) 

Since A^{H) is the sublist of aspeets of H prescribed by ^ and Ei^{E{H)) is 
the set of all sub-determined edges according to the sub-determination it fol¬ 
lows that {A^{H),E^{E{H))) is a MAG obtained from H aceording to the sub¬ 
determination As \A^{H)\ < \A{H)\, it follows that the order of M^(H) is 
lower than the order of H. Further, since self-loops may be created by edge sub¬ 
determination and discarded, and also sinee Ei^ is not injective, it follows that 
|£5(£(H))| < |£(H)|. 

2.3 MAG adjacency 

Two eomposite vertices are eonsidered adjacent if they share the same MAG edge, 
i.e. given two composite vertices u, v e V(//) are adjacent if and only if there is 
a MAG edge e e E{El) such that u, v e {Ko{e) ^Ttd{(^)] ■ Similarly, two MAG edges 
are considered adjacent if and only if they share a same eomposite vertex, i.e. two 
given edges e\,e 2 ^ E{H) are adjacent if and only if there is a composite vertex 
ue V(//) such that ue {Ko{e\),7^d{^i)} andue {^ 0 (^ 2 ), (^ 2 )}- 

Figureshows an illustrative example of three MAG edges. The figure depicts 
a four aspects MAG, where eaeh set of colored circles represents one aspect, and 
each edge has two elements of each aspect. 

The isolated edge (A1,A2,A3,A4,C1,C2,C3,C4) on the leftmost side of Fig¬ 
ure exemplifies the eomposite vertex adjaceney concept. In this ease the eom¬ 
posite vertices (A1,A2,A3,A4) and (C1,C2,C3,C4) are adjacent. The two edges 
(£l,£2,£3,£4,//l,//2,//3,//4) and {Hl,H2,m,HA,Kl,K2,K3,KA) exemplify 
a case of edge adjaeency. Sinee the composite vertex {Hl,H2, H3,HA) is shared 
by both edges, they are adjaeent. 

Although the strueture of a MAG edge is similar to an even uniform hyper¬ 
graph edge, the adjaeency definition used on MAGs is not the usual one adopted 
on hypergraphs. The adjacency concept found on a MAG is close to the one 
associated with traditional direeted graphs, where a MAG edge can be seen as a 
relation between two eomposite vertiees, whieh are eomposite objects constructed 
from aspect elements. Therefore, a MAG edge expresses a relationship between 
two (eomposite) objeets in the same way as a direeted graph edge. This concept 
leads to the isomorphism between MAGs and direeted graphs, as well as the elose 
relation between walks, trails, and paths on MAGs and directed graphs. 
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(A1,A2,A3,A4,C1,C2,C3,C4) 


(E 1,E2,E3,E4,H l,H2,H3,H4) (H 1,H2,H3,H4,K l,K2,K3,K4) 


Figure 1: Illustrative example of three MAG edges. 


2.4 MAG isomorphism 

In order to define MAG isomorphism it is neeessary to define the eoneept an aspeet 
list bijeetion. Given two MAGs H and K, both with p aspeets, an aspeet list bi- 
jeetion F :A{H) ~^A{K) is defined as a set of p bijeetive funetions, /i,/ 2 , 
sueh that eaeh aspeet of the MAG H is the domain of exaetly one of these fune¬ 
tions and eaeh aspeet of MAG K is the eodomain of exaetly one of these funetions. 
It follows from this definition that given a eomposite vertex u e Y{H), the aspeet 
list bijeetion F takes u to a eomposite vertex F’(u) e Y{K). 

Two MAGs of order p, H and K, are isomorphie if there is an aspeet list 
bijeetion F : A(H) A(K) sueh that an edge e e E(H) if and only if the edge 
{F{Me)),F{Me)))‘^E(K). 
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2.5 MAG walks, trails, and paths 


There is a close relation between walks, trails, and paths on a MAG and their 
counterparts in the isomorphic directed graph g{H). 

A walk on a MAG H is defined as an alternating sequence IT = [ui, ei, U 2 , £ 2 , U 3 , 

... .Uk-i^ek-i^Uk] of composite vertices u„ e Y{H) and edges ^ E{H), such 
that u„ = TToien) and u„+i = TZd{en) for 1 ^ n < k. It follows from this definition 
that in a walk, consecutive composite vertices as well as consecutive MAG edges 
are adjacent. 

We show in [19| that an alternating sequence IT of composite vertices and 
edges in a MAG // is a walk on H if and only if there is a corresponding walk 
Gw in the composite vertices representation of H. This means that a walk on a 
MAG H has a isomorphic walk on the directed graph g{H). Since trails and paths 
also are walks, we also show that the same isomorphism concept extends to them 
as well. 

Figure[^can also exemplify a MAG path. The two edges (£’1,£’2,£’3,£’4,//1, 
and m,HA,Kl,K2,K2,KA) can also be seen as part of 

the alternating sequence P = {E\,E2,E'i,EA), (£'l,£'2,£'3,£'4,//l,//2,//3,//4), 
{Hl,H2,m,m),{Hl,H2,m,HA,Kl,K2,K3,KA),{Kl,K2,K3,KA),w\)ic\\ char¬ 
acterizes a two-hops path from the composite vertex (El,£2,£3,£4) to the com¬ 
posite vertex (.^1,.^2,^3,^4). 

From the concept that walks, trails, and paths on a MAG have a isomorphism 
relation to their counterparts on the directed graph g{El), it follows that analysis 
and algorithms based on walks, trails, and paths can be formulated on the directed 
graph g{El). These properties will be extensively used in the current work. 


3 Algebraic Representation 


In this section, we discuss ways to represent MAGs [I9| by means of algebraic 
structures. As a consequence to the isomorphism between MAGs and traditional 
directed graphs, it is straightforward to construct matrix-based representations of 
MAGs. This section addresses these representations, using the MAG depicted in 
Figure as an illustrative example. 

Figurej^shows an example of a three aspect MAG T. It can be seen as the rep¬ 
resentation of a time-varying multilayer network, showing a small and simplified 
section of an urban transit system. More specifically. Figure [^depicts the MAG T 
in its composite vertices representation, g{T), which is the directed graph defined 
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in Expression ([^. 

Aligned with this view, the aspeets of MAG T ean be interpreted in the fol¬ 
lowing way: The first aspeet represents three distinet loeations, labeled 1, 2 and 3. 
Speeifieally, loeation 1 represents a subway station, loeation 2 a subway station 
with a bus stop, and loeation 3 a bus stop. The seeond aspeet represents two dis¬ 
tinct urban transit modes depicted as layers, namely Bus and Subway. Finally, the 
third aspect represents three time instants. The MAG edges can be seen with the 
following meaning: Location 1 has no edges on the bus transit mode, since it is a 
subway station. Similarly, location 3 has no edges on the subway mode, since it 
is a bus stop. The eight black edges represent bus and subway trips between loca¬ 
tions. As a simplification all trips are assumed to have the same duration. The red 
(dotted) edges represent the possibility of staying at a bus stop or subway station 
and not taking a transit. The six blue (dashed) edges show that it is possible to 
change between bus and subway layers at all times at location 2. As a simplifi¬ 
cation, the connection between the bus and subway layers is assumed to take no 
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time. We reeognize that the deeision of making these edges with 0 time length 
generates eyeles of length 0 in instanees of loeation 2. In real network analysis, 0 
length eyeles (and also negative length eyeles) ean eause problems. However, we 
ehoose to let these eyeles present in this toy example sinee they will eause no harm 
for the analysis eondueted in this thesis, and also, they make the toy example more 
eompaet and readable. Further, we remark that if desired, these 0 length eyeles 
eould be broken by adding new eomposite vertiees, or by making the subway/bus 
transition to have the same length as a subway/bus trip. 

In this model, walks represent the ways the urban transit system ean be used 
to travel from one loeation to another. For instanee, starting at loeation 1 on the 
subway layer at time tl, it is possible to reaeh loeation 3 on the bus layer at time 
t3. It ean be done by taking a subway trip to loeation 2 at time t2, switehing 
from subway to bus layer at loeation 2, time t2 and finally taking a bus trip from 
loeation 2 bus layer arriving at loeation 3 on the bus layer at time t3. 

The presenee of uneonneeted oeeurrenees of loeation 1 at bus layer and loea¬ 
tion 3 at subway layer ean be viewed as artefaets of the MAG eonstruetion. We 
eall these vertiees trivial eomponents of the MAG. This subjeet will be further 
addressed in this seetion. 

We remark that a Python implementation of all the algorithms presented in this 
seetion is available at the following URL: http : / /github . com/wehmuthklaus/ 
MAG_Algorithms, 

3.1 Companion tuple 

Although we show that every MAG H is isomorphie to a direeted graph designated 
g{H), it is important to note that the set of vertiees of this graph is Y{H), as shown 
in Expression (j^. Sinee the set Y{H) is the eartesian produet of all the aspeets 
in the MAG H, it is possible to reeonstruet the MAG’s aspeet list from Y{H), 
whieh is a step neeessary to obtain the MAG H from the direeted graph g{H). 

When the vertiees of the direeted graph G assoeiated with a given MAG H are 
not the eomposite vertiees themselves, it is neeessary to provide a meehanism to 
link eaeh vertex of the direeted graph to its eorresponding eomposite vertex on the 
MAG. This meehanism ean be, for instanee, a bijeetive funetion between Y{H) 
and y(G). 

In the eurrent work, we eonstruet representations for g{H), sueh as matriees, 
whieh do not direetly earry the tuples that eharaeterize the MAG’s eomposite ver¬ 
tiees. In this kind of representation, a vertex is assoeiated with a row or eolumn of 
a matrix. Therefore, additional information has to be provided in order to properly 
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link each row (column) of a matrix to its corresponding composite vertex on the 
MAG represented by this matrix. This is done by a bijective function D, defined 
in Section 3.2[ where D takes a composite vertex to a natural number, which is the 
row (column) number in the matrix. 

The implementation of D presented in this work is based on the concept 
of a companion tuple, which complements the matrix representation of a given 
MAG. For a MAG H with p aspects, its companion tuple has the form (|A(//)[1]|, 
\A{H)[2 ]\,..., \A{H)[p]\), so that the number of elements on it equals the order 
of H and each element represents the number of elements of an aspect of H. As a 
matter of notation, we represent the companion tuple of a given MAG H as 


l(H) = (|A(H)[l]|,|A(H)[2]|,...,|A(H)[p]|), 


(13) 


where p is the order of H. When there is no ambiguity in relation to which MAG 
we are referring to, we may use the notation T instead of t(//). For instance, the 
companion tuple of the MAG T shown in Figurej^is 'r(r) = (3,2,3), since T has 
3 aspects, of which the first has 3 elements, the second 2 elements, and the third 3 
elements. 

Algorithm shows the building of the companion tuple for a given MAG H. 
Assuming that the size of the aspect list A(//) and the size of each of the aspect 
sets contained mA{H) are known from the computational representation of A{H), 
the time complexity for building the companion tuple is 0 {p), where p is the 
number of aspects on MAG H. If, however, these sizes are unknown, then the 
time complexity is 0 {s), where 5 = Xif=i l^(^)[*]l’ since each element of each 
aspect has to be counted. We remark that, in either case, the time complexity for 
building the companion tuple is less than 0{Y{H)), which is the order of the set 
of composite vertices of the MAG. 


input :A{H) 


output: t{H) 


1 CompTuple(A(//)) 


2 

p^\A(,H)\ 

//number of aspects in the MAG 

3 

for i <— 1 to p do 


4 

n-i - WH)[>']i 

//number of elements in i-th aspect 

5 

end 


6 return T 



Algorithm 1: Construction of the companion tuple of a MAG. 
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For a given MAG H and a sub-determination we also define the sub¬ 
determined eompanion tuple whieh is obtained by multiplying eaeh entry 

of the original eompanion tuple by the equivalent entry of the tuple representation 
of as shown in Algorithm]^ The sub-determined eompanion tuple has the same 
value as the original eompanion tuple for the aspeets that have value 1 in ^ and 0 
otherwise. 


input 


output: 


1 SubCompTuple(T(//), 0 


2 

P- WH)I 

//number of aspects in the MAG 

3 

for i <— 1 to p do 


4 



5 

end 


6 return 



Algorithm 2: Construetion of sub-determined eompanion tuple. 


3.2 Order of composite vertices and aspects 


In general, the order of the eomposite vertiees and aspeets on a MAG is not rel¬ 
evant. That is, ehanging the order in whieh the aspeets or their elements are 
presented does not affeet the result of any algorithm or analysis performed on 
a MAG, sinee the MAG obtained by sueh ehanges is isomorphie to the original 
one. The definition of the MAG isomorphism adopted in this work ean be found 


in Seetion 2.4 However, in order to show the MAG’s algebraie representation in 
a eonsistent way, it is neeessary to link the MAG’s eomposite vertiees to rows and 
eolumns of matriees, whieh is aehieved by the bijeetive funetion D, defined in this 
seetion at Equation ( [T5| ). We now show the preliminary steps neeessary for the 
definition of funetion D, as implemented in this work. 

The aspeet order is adopted as the same in whieh the aspeets are plaeed on 
the MAG’s eompanion tuple. For the ordering of eomposite vertiees, we define 
the numerieal representation of eaeh eomposite vertex from its tuple. In order to 
obtain the eomposite vertex numerieal representation, we first translate the eom¬ 
posite vertex into a numerieal tuple. This is done by applying a family of indiees, 
one for eaeh aspeet on the eomposite vertex, where for every aspeet i the eorre- 
sponding index ranges from 0 to T; — 1, where T/ is the number of elements on the 
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i-th aspect of the MAG. Since this is a simple index substitution, we do not use a 
distinct notation for the composite vertex on its numerical tuple form. We, how¬ 
ever, reserve the notation u[z] to express the z-th element of the composite vertex 
on its numerical form. 

To calculate the numerical representation of a composite vertex, we define the 
weight of each position on the composite vertex tuple of a MAG H with p aspects 
as 


W(z,T) 


1 if z = 1, 

O/W Tj otherwise, 


(14) 


where z is the position in the tuple varying from 1 to p, T is the MAG’s compan¬ 
ion tuple, and Xj is the j-th element of the MAG’s companion tuple. Note that 
\x\= p meaning that the length of the companion tuple is the order of the MAG, 
i.e. the number of its aspects. Finally, we define the composite vertex numerical 
representation as 


D(u, x) = 1 + J^W(z, x) X u[z], 


(15) 


/=1 


where |t| = p is the MAG’s order, and v[z] is the z-th component of the com¬ 
posite vertex. Figure shows the MAG T with its composite vertices, and their 
numerical representations ranging from (1) to (18). In order to illustrate how the 
numerical representations are obtained, we show examples based on the MAG T. 
For this representation, we adopt aspect indices such that for aspect 1 we have 
Idx(l) = 0,Idx(2) = 1 and Idx(3) = 2. For aspect 2, we have Idx{Bus) = 0 and 
Idx{Subway) = 1, while for aspect 3, Idx{ti) = Q,Idx{t 2 ) = 1 and Idx{t'i) = 2. 
Since the companion tuple of MAG T is x{T) = (3,2,3) = x, the weights are 
VF(1 ,t) = 1, VF(2,t) = Ti = 3 and VF(3 ,t) = Ti x T 2 = 6. Therefore, the com¬ 
posite vertex v = {\,Bus,t\) has numerical representation D(v, T) = l-flxO-l- 
3x0-l-6x0 = 1, while D {{2, Subway, t 2 ), t) = 1-1-1x1-1-3x1-1-6x1 = 11 and 
D{{2,Bus,t2),x) = 1-1-1x1-1-3x0-1-6x2=14. 

Algorithm determines the numerical representation of a composite vertex 
V represented by its numerical tuple. The presented implementation extends the 
concepts presented in Equations ( fK] ) and ( [T5j ), so that this algorithm can also 
be used to determine the numerical representation of sub-determined composite 
vertices. In order to determine the numerical representation of a sub-determined 
vertex, function D shown in Algorithm receives the full composite vertex tu¬ 
ple (not sub-determined) and the sub-determined companion tuple. The if seen 
at line 7 of Algorithm makes that the 0 entries found in a sub-determined com¬ 
panion tuple are discarded for the construction of the sub-determined numerical 
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representation of the composite vertex. The time complexity for this algorithm is 
0{p), where p is the number of aspects on the MAG in question. 

Given the numerical representation of any composite vertex, it is possible to 
reconstruct its tuple. In order to do this, we calculate the numerical value of the 
index of each element on the tuple, as 

N{dp,x) = [{{d-\)modW{i+\,x)) /1T(z,t)J, (16) 

where d is the composite numerical representation, i is the position of the com¬ 
posite vertex tuple to be calculated, x is the MAG’s companion tuple, mod is 
the modulus (division remainder) operation and [xj is the floor operator, which 
for any x e M corresponds to the largest integer i e Z such that i < x. Note that 
for calculating N{d,p, x) for a MAG with p aspects, it is necessary to calculate 
W{p + l,x). Considering the definition of W from Equation ( [T^ , it follows that 
W{p+l,x) = 0^=1 = l^(^)l’ the number of composite vertices on the MAG. 

For instance, taking the composite vertex with numerical representation 14 of 
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input : V, x{H) 
output: d 

1 

2 D(v,'r(//)) 


P 

d 


m)\ 

0 

W <— 1 

for / <— 1 to do 

if T (//)[/] 7 ^ 0 then 
d d + (v[/] * w) 
w <— w* t (//)[/] 

end 
end 


3 

4 

5 

6 

7 

8 
9 

10 
11 

12 return d 


// V is the numerical tuple of the composite vertex 
//number of aspects in the MAG 


Algorithm 3: Determination of the numerieal representation of a eom- 
posite vertex. 


the MAG T, we have that 

A(14,l, (3,2,3)) =[((14-l)moJ3)/lJ = [1/lJ = 1 
A(14,2, (3,2,3)) =[((14-l)moJ 6)/3j = [1/3J = 0 
A(14,3, (3,2,3)) =[((14-l)moJ 18)/6J = [13/6J = 2. 


We ean therefore define the inverse of funetion D as 

D-\d, x) = iN{d, 1, t),A(J,2, t), ... ,A(J, |t|, t)), (17) 

whieh reeonstruets the eomposite vertex tuple in its numerieal form. From this, 
we ean see that, for instanee, D“^(14, (3,2,^) = (1,0,2), whieh eorresponds to 
the eomposite vertex {l^Busf'i). Algorithm^ shows the implementation of . 


The relation between the eomposite vertex numerieal representation and its 
tuple ean also be seen as a eonsequenee of the natural isomorphism between the 
MAG H and its eomposite vertiees representation, g{H). The role of this relation 
will beeome elear in Seetions [T4| to where the matrix forms of the MAG are 
presented. 
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input : d, x{H) 
output: V 

1 

2 InvD((i, x{H)) 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 return v 


P^\r(H)\ 
w[l] <— 1 
w/ <— 1 

for / <— 1 to do 

if x{H)[i] ^ 0 then 

w[i + 1] <— w\i\ * T{H)[i] 
w/ <— wZ + 1 

end 
end 

for Z <— 1 to wl do 

I v[Z] <— {d Mod w\i + l])/w[Z] 

end 


// V is the numerical tuple of the composite vertex 
//number of aspects in the MAG 


Algorithm 4: Determination of the eomposite vertex from its numerieal 
representation. 


3.3 Elimination of trivial components 

In the MAG T shown in Figure the eomposite vertiees of numerieal represen¬ 
tation (1), (6), (7), (12), (13), and (18) are trivial eomponents (i.e. uneonneeted 
eomposite vertiees). They are ereated in eonsequenee of the regularity needed on 
the MAG H to build the set N{H). This type of padding is not neeessary in a 
direeted graph and its algebraic representation. Therefore, it is possible to remove 
the trivial components from the composite vertices representation and its associ¬ 
ated matrices. However, it is important to bear in mind that the graph resulting 
from this transformation may no longer be isomorphic to the MAG and neither 
are the matrices associated with it. The only case in which the isomorphism is 
preserved is when there are no trivial components on the MAG and nothing is 
removed. Nevertheless, this kind of transformation can be helpful for applica¬ 
tion, by reducing the number of composite vertices present on the graph and so 
simplifying its construction and manipulation. 


The same sort of padding is discussed in [ 16|, where authors suggest that this 
padding may cause problems in the computing of some metrics, such as mean 
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degree or elustering eoeffieients, unless one aeeounts for the padding seheme in 
an appropriate way. In this subseetion, we show that the padding with the triv¬ 
ial eomponents may be eliminated, if desired. Anyway, if needed, it suffiees to 
be eautious in eomputing the metries of interest on MAGs by eonsidering the 
existenee of the padding seheme, as suggested by p^ . In partieular, the MAG 
algorithms we diseuss in Seetion]^ remain unaffeeted by this padding issue. 

For a given MAG H, we define its main eomponents graph m{H) as the MAG’s 
eomposite vertiees representation with all its trivial eomponents removed. Fig- 
ure|^shows the main eomponents graph m{T) for the MAG T. It is worth noting 
that numerieal representations are not defined for m(r). 



This ean be aehieved algebraieally for any MAG H with the help of a matrix 
R(//) eonstrueted from the identity I e where n = |V(//)| is the number 

of eomposite vertiees on the MAG. The matrix R(//) is obtained from this nxn 
identity by removing the eolumns whieh mateh the numerieal representations of 
the trivial eomponents of the MAG. Therefore, assuming that the MAG H has r 
trivial eomponents, the matrix R(//) e ^ rows and n — r eolumns. In 

partieular, in the eases where the MAG H has no trivial eomponents, we have that 
R(//) 

It is also worth noting that the matrix Im(^) = R(^) R(^) 
trix akin to the identity I e but the diagonal entries eorresponding to the 

trivial eomponents (removed in R(//)) have value 0. Therefore, multiplying a 
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ny. n matrix by !„(//) to the left has the effeet of turning all entries on the rows 
eorresponding to the trivial eomponents to Os. Similarly, multiplying by \m{H) to 
the right has the effeet of turning the entries of the eolumns eorresponding to the 
trivial elements to Os. 

As an example, we show the matrix R(r) e 


R(r) 


0 0 
1 0 
0 1 

0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
1 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
1 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
1 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
1 0 
0 1 
0 0 
0 0 
0 0 


0 o' 

0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 ’ 
0 0 
0 0 
0 0 
0 0 
0 0 
1 0 
0 1 
0 0 


( 18 ) 


whieh is obtained from the 18 x 18 identity matrix by removing the eolumns 
1,6,7,12,13, and 18 that eorrespond to the trivial eomponents of the MAG T. 


3.4 Adjacency matrix 


As a direet eonsequenee of the isomorphism between MAGs and traditional di- 
reeted graphs, it is expeeted that a MAG ean be represented in matrix form. In 
faet, sueh representations ean be aehieved direetly by the eomposite vertiees rep¬ 
resentation of MAGs, presented in Seetion 2.1[ Sinee for any given MAG H its 
eomposite vertiees representation is a traditional direeted graph, it ean be repre¬ 
sented in matrix form. 

One of sueh representations is the MAG’s adjaeeney matrix. This matrix is 
obtained from the MAG’s eomposite vertiees representation, g{H), and its eom- 
panion tuple t(//). In faet, the MAG’s adjaeeney matrix is the adjaeeney matrix 
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of the composite vertices representation, where the order of the rows and columns 
is given by the numerical representation of the composite vertices of g{H). 

Since the set Y{H) of composite vertices of a given MAG H is obtained by 
the cartesian product of all aspects of the MAG (as shown in Expression Q), it 
follows that the number of composite vertices on a given MAG H with p aspects 
is 

p 

„=|V(//)| = f]T,-, (19) 

(=1 

where T, is the z-th element of the MAG’s companion tuple, i.e. the number of 
elements on the MAG’s z-th aspect. 

The general form of any entry of the matrix J {H) is given by 


I 1 if(u,v)eE(g(//)), 
[ 0 otherwise, 


( 20 ) 


where (u,v) e E{g{H)) means that (u,v) is an edge on the composite vertices 
representation g{H) of the MAG H, so that u,v e Y{H) are composite vertices 
of H. It follows from the definition of g{H) and its natural isomorphism to H, 
that (u,v) e E{g{H)) if and only if there is an edge e e E{H) such that u = TToie) 
and v = It is important to note, however, that the notation y’u.v is in fact 

a shorthand for 7£)(u,t),d(v,t)’ where T>(u, r) is the row number and T>(v, r) the 
column number of the matrix entry. This ties the construction of the adjacency 
matrix of a MAG with its companion tuple, since it is used in the determination 
of the numerical representation of a composite vertex (D(u, t)). Therefore, the 
adjacency matrix of any given MAG is always presented with its companion tuple. 

The adjacency matrix of a given MAG H is constructed by Algorithmj^ where 
|V(//)| is the number of composite vertices in H, which can be calculated using 
Equation ( [T^ , D(;ro(e), t) and D(;r^(e), t) are the numerical representation of the 
origin and destination composite vertices of edge e e E{H), respectively, as de¬ 
fined in SectionConsidering that a sparse matrix with all entries 0 can be cre¬ 
ated in constant time, and that both functions CompTuple and D (see Algorithm[^ 
and Algorithm]^ have time complexity 0{p),we conclude that Algorithm]^ has 
time complexity 0{p * \E where p is the number of aspects of MAG H and 
\E{H)\ the number of edges. 

As an example, the adjacency matrix of the MAG T is shown in Expres¬ 


sion (21). This adjacency matrix J(r) 
22 are non-zero. 


>18x18 


has 324 entries, of which just 
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input :H={A,E) 
output: 

1 AdjMatrix(//) 

2 n^|V(//)| 

3 T CompTuple{A{H)) 

4 J(H) n X n matrix with all entries = 0 

5 for each e e E (H) do 

6 D{noie),T) 

7 v<~D{7td{e),T) 

8 J(iy)[u,v] <- 1 

9 end 
ro return 


//companion tuple 


// numerical origin 
// numerical destination 


Algorithm 5: Building J{H) from MAG H. 


J(7’) 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 10 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
0 10 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
0 1 1 
0 1 1 

0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 10 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
0 0 0 
0 0 0 
1 1 0 
1 1 0 
0 0 0 
0 0 0 
0 10 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 1 1 
0 1 1 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 10 
0 0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
1 1 0 
1 1 0 
0 0 0 
0 0 0 
0 10 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


(21) 


It is important to note that the order of the columns and rows of J(7’) is 
given by the numerical representation of the composite vertices. Thus, for in¬ 
stance, the 1 at row 2, column 8 represents the edge between the composite ver¬ 
tices with numerical representations 2 and 8, witch in turn represents the edge 
{2,Bus,t\,2,Bus,t2) of the MAG T. In this way, although J(r) is presented in 
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matrix form, together with the eompanion tuple T(r), it fully represents the MAG 
T, earrying the proper adjaeeney notion used to define transitive eonstruetions, 
sueh as walks and paths on the MAG. 

For an arbitrary MAG H, its main eomponents graph m{H) is obtained by 
removing the MAG’s trivial eomponents, as stated in Seetion 3.3[ The matrix 
J(m(//)) is then obtained with the use of the matrix R(//), presented in See¬ 
tion 


33| is obtained as 

J(m(//)) =R(//)^ J(//) R(//), 


( 22 ) 


where J(m(//)) e ^n-r^n-r adjaeeney matrix eontaining only the main eom¬ 
ponents of the MAG. 

It is also possible to obtain the adjaeeney matrix i{H) from J(m(//)). This 
follows from the faet that on the adjaeeney matrix i{H) the rows and eolumns 
eorresponding to trivial eomponents are already zero. Therefore, 


J(//) =!,„(//) (23) 

where \m{H) = R(//) R(//)^. Then, we have that 

R(//) J(m(//)) K{hY = R(//) i{H) K{H) (24) 

= \^{H) J{H) !,„(//) 


Expression (25) shows J(m(r)), the adjaeeney matrix of m{T). This matrix 
is obtained from the adjaeeney matrix J(r) by removing the rows and eolumns 
whieh represent the trivial eomponents of the MAG T. In this ease, the trivial eom¬ 
ponents are the eomposite vertiees with numerieal representations 1,6,7,12,13, 
and 18. This matrix is ealeulated as J(m(r)) = R(r)^ J(r) R(r), so that 


J(m(r)) 


0 0 
0 0 
0 0 
1 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 1 

0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


1 1 
1 1 

0 0 
0 0 
0 0 
0 0 
0 0 
1 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
1 1 
1 1 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
1 1 
1 1 
0 0 
0 0 
0 0 
0 0 
0 0 
1 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
1 1 
1 1 
0 1 
0 0 
0 0 
0 0 


(25) 
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In general, the adjaeeney matriees assoeiated with a MAG are sparse, meaning 
that for an R X n adjaeeney matrix the number of non-zero entries of the matrix 
is of the order 0{n). Sinee the non-zero entries on the MAG adjaeeney matrix 
eorresponds to the edges present on the MAG, the adjaeeney matrix being sparse 
means that the number of edges m on the MAG is of the same order of the number 
of eomposite vertiees n, i.e. m is of order 0{n). Therefore, these matriees ean 
be stored effieiently using sparse matriees representations, sueh as Compressed 


Sparse Column (CSC) or Compressed Sparse Row (CSR) p4| . Assuming that 
the number of edges is larger than the number of eomposite vertiees, these repre¬ 
sentations provide a spaee eomplexity of 0{m) for storing the MAG’s adjaeeney 
matriees. Further, they also provide effieient matrix operations, whieh will be 
explored in the algorithms presented in Seetion]^ 


3.5 Incidence matrix 

Given that every MAG is isomorphie to a direeted graph, it follows that it ean 
be represented by an ineidenee matrix (and its eompanion tuple). For any given 
MAG H, this matrix is eonstrueted from the eomposite vertiees g{H) and the 
eompanion tuple x{H), adopting the vertex order indueed by the numerieal repre¬ 
sentation presented in Seetion 3.2 The MAG’s ineidenee matrix C(//) e 
where m = \E{H)\ is the number of edges in the MAG and n = |V(//)| is the 
number of eomposite vertiees on the MAG, is defined then as 


^e,u ~ 


1 if\x = no{e), 
-1 \f\x = Kd{e), 
0 otherwise, 


(26) 


where e e E{g{H)) is an edge in MAG H and u e Y{H) is a eomposite vertex 
in MAG H. Here, the notation Ce,u is a shorthand for c/^(e),D(u,T)> where Id{e) 
is an numerieal index for eaeh edge and D{u, t) is the numerieal representation 
of the eomposite vertex u. Note that the use of the eomposite vertex numerieal 
representation ties the ineidenee matrix to the MAG’s eompanion tuple. 

Although the order of the eomposite vertiees is defined by eaeh eomposite 
vertex numerieal representation, the order used to represent the MAG edges in 
the ineidenee matrix is not relevant. The ineidenee matrix of a direeted graph has 


several well-known properties [351, among whieh, the property that the ineidenee 
matrix of a direeted graph with k eonneeted eomponents has rank n — k, where n 
is the number of vertiees of the graph. This property is useful for defining other 
matriees based on the ineidenee matrix. 
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For a given MAG H, the ineidence matrix C(//) is built by Algorithm]^ where 
D{Ko{e),T) and D{nd{e),T) are the numerieal representation of the origin and 
destination eomposite vertiees of edge e e E{H), respeetively, as defined in See- 
tion 3.2 and Id{e) is a unique numerieal index for the edge e e E{H), ranging 


from 1 to m. Considering that a sparse matrix with all entries 0 ean be ereated in 
eonstant time, and that both funetions CompTuple and D (see Algorithm [T] and 
Algorithm]^ have time eomplexity of 0{p), we eonelude that Algorithm]^ has 
time eomplexity of 0{p* \E{H)\), where p is the number of aspeets of MAG H 
and \E{H)\ the number of edges. 


input :H={A,E) 
output: C{H),t:{H) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


IncidMatrix(//) 

|£(H) 

T CompTuple{A{H)) 

C{H) <— mx n matrix with all entries 
for each ee E{H) do 

u<- D{Ko{e),T) 
y ^ D{nd{e)J) 

C{H)[i,u] <- 1 
C{H)[i,y]^-l 


12 


end 


13 return C(//),r 


0 


//companion tuple ofH 


// index of edge e 
// numerical origin 
// numerical destination 


Algorithm 6: Building C{H) from MAG H. 


Given the incidence matrix C{H) of a MAG H, it is possible to obtain the 
incidence matrix of the main components graph (Z{m(H)) using the matrix R(//) e 
defined in Section 3.4 The incidence matrix of m{H) is given by 


C(m(//)) = C{H) K{H) 


(27) 


Further, given the incidence matrix of the MAG’s main components graph and 
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the matrix R(//), it is possible to reeover the MAG’s ineidenee matrix, as 

C(m(//)) R{Hf = C{H) R{H) R{Hf (28) 

= C{H) l^H) 

= C{H). 

This is only possible beeause the eolumns of C{H), whieh are foreed to 0 by the 
multiplieation by were already 0, as the eomposite vertiees represented by 

them have no edges ineident to them. 

The ineidenee matrix C(r) of the example MAG T is shown in Expression (j2^. 
The vertiees (eolumns) order is determined by the vertiees numerieal representa¬ 
tion, while the edge order remains uneonstrained. The trivial eomponents eorre- 
spond to eolumns 1,6,7,12,13, and 18, whieh have all entries with value 0. 


C{T) 


0 1 0 

0-10 
0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 10 

0 0 1 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 1 0 

0 0 1 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


0-10 
0 1 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

1 0 0 

0 1 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

1 0 0 

0 1 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


0 0 0 

0 0 0 

0 10 

0-1 0 

0 0 0 

0 0 0 

0-1 0 

0 0-1 

0 0 0 

0 0 0 

0 10 

0 0 1 

0 0 0 

0 0 0 

0 0-1 

0-1 0 

0 0 0 

0 0 0 

0 10 

0 0 1 

0 0 0 

0 0 0 


0 0 0 

0 0 0 

0-10 
0 1 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

-1 0 0 

0-10 
0 0 0 

0 0 0 

1 0 0 

0 1 0 

0 0 0 

0 0 0 

0-10 
-1 0 0 

0 0 0 

0 0 0 

1 0 0 

0 1 0 


0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 10 

0-1 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0-1 0 

0 0-1 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0-1 

0-1 0 

0 0 0 

0 0 0 


0 0 0 - 

0 0 0 

0 0 0 

0 0 0 

0-10 
0 1 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

-1 0 0 

0-10 
0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0-10 
-1 0 0 _ 


(29) 


The main eomponents ineidenee matrix C(m(r)) is depieted in Expression ( [30] ). 
It is obtained from the matrix C(r) by removing the trivial eomponents, as shown 
in Expression ([28|). 
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C(m(r)) 


100 -1 00000000 
-1 00100000000 
0000100 -1 0000 

0000 -1 0010000 

00000000100 -1 

00000000 -1 001 

1000 -1 0000000 

01000 -1 000000 

001000 -1 00000 

0001000 -1 0000 

00001000 -1 000 

000 0 0 1 0 0 0 -1 0 0 

000 0 0 0 1 0 0 0 -1 0 

00000001000 -1 

10000 -1 000000 

0100 -1 0000000 

0010000 -1 0000 

000100 -1 00000 

000010000 -1 00 
00000100 -1 000 

00000010000 -1 

0000000100 -1 0 


(30) 


In general, the incidence matrices related to MAGs are sparse, and therefore 
can be efficiently stored using sparse matrices representations, such as CSC or 
CSR [ (^ . Assuming that the number of edges on the MAG is larger than the 
number of composite vertices, the use of these representation lead to a memory 
complexity of 0{m), where m= \E(H) \ is the number of edges on the MAG H. 

3.6 Laplacian matrices 

3.6.1 Combinational Laplacian 

We construct the combinational Laplacian matrix of a given MAG H from its 
incidence matrix C{H), as 



(31) 


Since C{H) is an m x n matrix, it follows from this construction that, as expected, 
the Laplacian L(//) is a positive semidefinite nx n matrix. Further, since the rank 
of C{H) is n — k, where k is the number of connected components of H, it follows 
that the rank of L(//) is also n — k. Consequently, the dimension of the nullspace 
of L(//) is k, the number of connected components on the MAG H, a well-known 
property of the Laplacian matrix. 

In the case of the Laplacian L(//), each one of the trivial components of the 
MAG counts as a distinct connected component. Therefore, for a MAG with t 
trivial components, we have that k^t, the equality happening in the case where 
the MAG only has trivial components, i.e. when the MAG has no edges. 
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The Laplacian L(r) of the example MAG T is given by 


L(r) = 


0 

4 

0 

0 

-2 

0 

0 

-1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

-2 
0 

0 0 0 

4 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

-10 0 

-10 0 

0 0 0 

0 0 0 

0 
0 
0 
0 
0 


0 0 0 

0 0-1 

0 0-1 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 6 

0 0 0 

0 

0 0-2 

0 0 0 

0 0 0 

0 0-1 

0 0-1 

0 0 0 

0 0 0 

0 0 0 


0 

0 

0 

-1 

-1 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

-1 

-1 

0 


0 

0 

0 

-1 

-1 

0 

0 

-2 

0 

0 

6 

0 

0 

0 

0 

-1 

-1 

0 


0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 - 1 - 

0 0 - 1 - 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 4 

0 0 0 

0 0 0 

0 0-2 

0 0 0 


(32) 

Sinee the sum of all eolumns of L(r) is 0 and six of the eolumns are 0, it 
follows that the dimension of the nullspaee of L(r) is 7, whieh is the expeeted 
value, as the MAG T has 6 trivial eomponents and a single main eomponent. The 
entries with value —2 refieet the faet that in this direeted graph there are pairs of 
opposing direeted edges, whieh ean be interpreted as a bi-direetional eonneetion. 


As is shown in Seetion 3.6.2 this ean also be seen as the weight assoeiated with 
this eonneetion. 

The Laplaeian ean also be eonstrueted for the main eomponents of a given 
MAG H. In this ease, the Laplaeian is eonstrueted as 


or 


L(m(//)) = K{Hf L(//) R(//), 
L(m(//)) = 

The main eomponent Laplaeian for the MAG T is 


(33) 


(34) 


L(m(r)) = 


4 

0 

0 

-2 

-1 

-1 

0 

0 

0 

0 

0 

0 


-2 

0 

0 

4 

0 

0 

-1 

1 


-1 

-1 

0 

0 

6 

0 

0 

-2 


0 -1 - 


0 

0 

0 


1 

0 

0 


0 

0 

-1 

-1 

-2 

0 

0 

6 

0 

0 

-1 

-1 


0 

0 

0 

0 

-1 

-1 

0 

0 

4 

0 

0 

-2 


0 

0 

0 

0 

0 

0 

-1 

1 


0 -2 

0 0 


0 

4 


(35) 


Sinee the six trivial eomponents were eliminated, the dimension of the nullspaee 
of L(m(r)) is 1. 
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3.6.2 Weighted Laplacian 


The weighted Laplaeian matrix of a MAG H is obtained in a similar way to the 
eombinational Laplaeian. However, an additional diagonal weights matrix is used 
to assoeiate a weight to eaeh of the edges of H. We denote a weights matrix for a 
given MAG // as W (H) e where m = \E(H) \ is the number of edges in H. 

Given a MAG H and a weights matrix W{H), the weighted Laplaeian is defined 


as 


= C{H). 


(36) 


In general, the entries on the main diagonal of a weights matrix W(//) are pos¬ 
itive real values. In this ease, W(//) is a symmetrie positive-definite matrix and, 
therefore, .^{H) is symmetrie positive-semidefinite. Henee, the rank of .^{H) is 
the same as the rank of C(//), so that the nullspaee of has the same dimen¬ 
sion as the nullspaee of L(//). It ean be seen that this matrix represents the same 


objeet as the supra-laplaeian deseribed in [36|. Nevertheless, here it is obtained 
direetly from the MAG’s representation by matrices and further, distinct weights 
can be directly assigned to each edge if the application needs it. 

As an example of weighted Laplacian for the MAG T, consider a weight ma¬ 
trix W(r), where the values of the entries on the main diagonal are given by 


Diag(W(r)) = [0.5,0.5,0.5,0.5,0.5,0.5,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]. 


(37) 

This weights matrix assigns weight 0.5 to all six edges that form the bidirectional 
connection between layers on the example MAG T. This effectively converts this 
edge pairs into a undirected edge. By doing this, the obtained weighted Laplacian 
matrix has the more familiar structure associated with the Laplacian of undirected 
graphs. For this weights matrix, we have 


^{T) 


0 0 

0 3 

0 0 

0 0 

0 -1 

0 0 

0 0 

0 -1 

0 -1 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 


0 0 

0 0 

2 0 

0 2 

0 0 

0 0 

0 0 

-1 0 

-1 0 

0 -1 
0 -1 
0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 


0 0 

-1 0 

0 0 

0 0 

3 0 

0 0 

0 0 

0 0 

0 0 

~1 0 

-1 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 


0 0 

0 -1 

0 -1 

0 0 

0 0 

0 0 

0 0 

0 5 

0 0 

0 0 

0 -1 

0 0 

0 0 

0 -1 

0 -1 

0 0 

0 0 

0 0 


0 0 0 0 

-10 0 0 

-10 0 0 

0-1-10 
0-1-10 
0 0 0 0 

0 0 0 0 

0 0-10 
4 0 0 0 

0 4 0 0 

0 0 5 0 

0 0 0 0 

0 0 0 0 

-10 0 0 

-10 0 0 

0-1-10 
0-1-10 
0 0 0 0 


0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 - 1-1 0 

0 - 1-1 0 

0 0 0 -1 

0 0 0 -1 

0 0 0 0 

0 0 0 0 

0 3 0 0 

0 0 2 0 

0 0 0 2 

0-100 
0 0 0 0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

-1 

0 

0 

3 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


( 38 ) 
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If considering only the main eomponents m{H) of a given MAG H, we have 
= R{Hf C{Hf W{H) C{H) R{H) = R{Hf ^{H) R{H). (39) 


For the ease of the example MAG T and the weights matrix deseribed by 
Equation ([37]), 


^(m(r)) 


3 0 

0 2 

0 0 

-1 0 

-1 -1 

-1 -1 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 


0 -1 

0 0 

2 0 

0 3 

0 0 

0 0 

-1 -1 

-1 -1 

0 0 

0 0 

0 0 

0 0 


-1 -1 

-1 -1 

0 0 

0 0 

5 0 

0 4 

0 0 

-1 0 

-1 -1 

-1 -1 

0 0 

0 0 


0 0 

0 0 

-1 -1 

-1 -1 

0 -1 

0 0 

4 0 

0 5 

0 0 

0 0 

-1 -1 

-1 -1 


0 0 

0 0 

0 0 

0 0 

-1 -1 

-1 -1 

0 0 

0 0 

3 0 

0 2 

0 0 

-1 0 


0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

-1 -1 

-1 -1 

0 -1 

0 0 

2 0 

0 3 


(40) 


3.6.3 Normalized Laplacian 


Another form of applying weights to the Laplaeian matrix on a given MAG H 
leads to the equivalent of the normalized Laplaeian matrix [37|. In this ease, 
the weights are applied to the eomposite vertiees instead of the edges, as in the 
weighted Laplaeian. In order to obtain the normalized Laplaeian, weights are 
applied to the non-zero eolumns of the ineidenee matrix C{H), whieh eorrespond 
to the eomposite vertiees of H that are not trivial eomponents (i.e. uneonneeted 
eomposite vertiees). The weights applied to the non-zero eolumns are sueh that 
the veetor represented by eaeh eolumn beeomes an unitary veetor. This leads to a 
diagonal weights matrix N{H) e where n = |V(//)|, for whieh 


N(H)/j - 


l/||c,|| if / = 7, and eolumn i A 0 , 
0 otherwise, 


(41) 


where c, is the z-th eolumn of the ineidenee matrix C{H) and ||c/|| is the Euelidean 
norm of Ci. The normalized Laplaeian is then obtained by 

£(//) = {C{H) {C{H) N(//)) (42) 

= N(//)^ C{H) N(//) 

= N(//) L{H) N(//). 
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Since ||c/|| = Ij-s/cfi, where di is the degree of the eomposite vertex eorresponding 
to eolumn Ci, it follows that the formulation for shown in Equation (42) 
eoineides with the one proposed in p7|. 

As with the other kinds of Laplaeian matriees, the trivial eomponents of 
can be eliminated using the matrix R(//) as 


£.{m{H)) = R{Hf £(//) R{H) 


(43) 


4 MAG Algorithms 


The MAG algorithms eovered in this seetion are based on the MAG’s adjaeeney 
matrix or on its adjaeeney list. Sinee in general we expeet the adjaeeney matrix 


to be represented using sparse CSR ou CSC formats p4| , it follows, due to the 
strueture of the CSR and CSC formats, that the adjaeeney matrix and adjaeeney 
list ean be seen as very elosely related representations. The algorithms used in 
MAGs are direetly derived from the basie well-known algorithms used with di- 
reeted graphs In this sense, the purpose of this seetion is not to 

propose new algorithms, but to show how known algorithms may be adapted for 
applieation in MAGs. 

We remark that a Python implementation of all the algorithms presented in this 
seetion is available at the following URL: http : / /github . com/wehmuthklaus/ 
MAG_Algorithms, 


4.1 Auxiliary matrices and vectors 

When operating upon a matrix representation, a few auxiliary matriees and veetors 
are neeessary to express the desired operations. We now define these veetors, 
whieh are used on the remainder of this seetion: 

1. All Os 

We denote 0 the eolumn veetor with all entries equal to 0. Usually we as¬ 
sume that 0 has the right dimension (i.e. number of rows) for the indieated 
operation. When neeessary to improve readability, we indieate the dimen¬ 
sion by sub-seript as in 0„. 

2. All Is 

We denote 1 the eolumn veetor with all entries equal to 1. Usually we as¬ 
sume that 1 has the right dimension (i.e. number of rows) for the indieated 
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operation. When neeessary to improve readability, we indieate the dimen¬ 
sion by sub-seript as in 1„. 

In all oases, we assume the veotors have the dimension neeessary for the oper¬ 
ation where they is applied. 

Moreover, speeially oonstruoted matrioes are used to build sub-determined 
algebraio algorithms for MAGs. These matriees provide reduotion/aggregation 
operations needed for sub-determined algorithms. Although these matriees are 
speeially oonstruoted for the MAG and the sub-determination in question, they 
have distinct properties and can be constructed by a general algorithm. In fact, the 
construction of sub-determined algorithms relies on the use of functions to aggre¬ 
gate/reduce results according to the applied sub-determination. In some cases, this 
function can be as simple as just summing up values obtained in composite ver¬ 
tices, which are reduced to the same sub-determined vertex. However, depending 
on the algorithm being constructed, this aggregation may need a more elaborate 
function, which may not be expressed in terms of matrix multiplications. 

Given a MAG H and a sub-determination the sub-determination matrix 
M|r(//) e is a rectangular matrix, where n = |V(//)| is the number of com¬ 
posite vertices of H and m = |V^(//)| is the number of composite vertices of 
the sub-determination ^ applied to the MAG H. Since a sub-determination is a 
(proper) subset of the aspects of a MAG, it follows that m\n, i.e. the number of 
composite vertices of a MAG is a multiple of the number of composite vertices in 
any of its sub-determinations. Further, has the property of having exactly 

one non-zero entry in each column, and the position of this entry is determined by 
the numerical value of the sub-determined composite vertex. 

Algorithm [tI shows the construction of the sub-determination matrix (//) 
for a given MAG H and sub-determination The function D takes a compos¬ 
ite vertex to its numerical representation and the function 5^ takes a composite 
vertex to its sub-determined form, i.e. it drops the aspects not present in the sub¬ 
determination. To determine the time complexity of Algorithm |7} we consider 
that the count of composite vertices in line 3 is 0(|V(//)|), the same is the case 
for the count on line 4, the construction of companion tuple at line 2 is 0{p), the 
construction of an empty sparse matrix at line 5 is 0(1), and, finally, the for loop 
initiated at line 6 is also 0(|V(//)|). Since the number of aspects p « |V(//)|, we 
conclude that the time complexity of Algorithm|^is 0(|V(//)|). 

For instance, consider the example MAG T and a sub-determination = 
OII 2 , which drops the third aspect of T. The aspect dropped is the aspect of 
time instants and, therefore, the two aspects present in Q are location and tran- 
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input : T (H) and ^ 

output: M^{H) 

1 SubDetMatrix(T(//), Q 

= SubCompTuple{x{H), Q 

»-|V(H)| 

m <— |V^(//)| 

<— mx n sparse matrix 
for 7 ■«— 1 to n do 

u^D-1(7,t(//)) 
i^D{u,T^) 

end 

11 return (H) 


2 

3 

4 

5 

6 

7 

8 
9 

10 


A sub-determined companion tuple 


//numeric tuple form ofj 
// sub-determined numerical representation 


Algorithm 7: Construction of M^. 


sit layers. Since in T there are 3 locations and 2 transit layers, it follows that 
|V^f(r)| = 6. Hence, M|rj(r) e constructed according to Algorithmis 

given by 


Mc,(r) 


1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 


As a further example, consider the MAG T and a sub-determination = 
IOO 2 , which drops the location and transit layer aspects, leaving only the time 
instants aspects. Since there are 3 time instants in T, it follows that e 

^3x18 js 


M^AT) 


111111000000000000 

000000111111000000 

000000000000111111 


(45) 


Note that in these cases the multiplication by the sub-determination matrices per¬ 
forms the sum of the distinct composite vertices that are reduced to a same sub¬ 
determined vertex. For instance, given the sub-determination = 01 12 , the ma¬ 
trix M^,(r) is used to aggregate values found in 3 composite vertices into a single 
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sub-determined vertex. The aggregation funetion in this ease is a simple sum. 
The same is done by the matrix M|ry(r) for the sub-determination = IOO 2 , 
where in this ease eaeh sub-determined result is the sum of values obtained for 6 
eomposite vertiees. 


4.2 Universality of matrix algorithms 


In this seetion, we show that every funetion that ean be obtained from a MAG 
to a given eo-domain set ean also be obtained from a matrix representation of 
the MAG. Here the set El is the quotient set of finite MAGs under isomorphism 
defined in Seetion |2.1[ Note that a permutation cr of a given adjaeeney matrix 


i{H), together with the funetion D(j, represents the same MAG H as i{H), so 
that permutations of adjaeeney matriees are isomorphie. Thus, we have the set 
J, whieh is a quotient set of pairs of adjaeeney matriees and assoeiation 

funetions D, under adjaeeney matrix permutations. Therefore, an element of J is 
an equivalenee elass of adjaeeney matriees and D funetions. Sinee we eonsider the 
pair t(//)) as the eanonieal adjaeeney matrix representation of the MAG 

H, we assign this pair as the elass representative of the MAG H in J. 


Theorem 1. The adjacency matrix J {H) and companion tuple t{H) obtained from 
the MAG H by Algorithm^are isomorphic to the MAG H. 

Proof. We show that Algorithm ean be seen as a funetion that takes a given 
MAG H to its adjaeeney matrix and eompanion tuple, and that this funetion pre¬ 
serves the adjaeeney strueture of the original MAG. Further, we show that, from 
the adjaeeney matrix i{H) and eompanion tuple x{H), we ean eonstruet a MAG 
H that is isomorphie to MAG H. 


Given the sets El and J, Algorithm]^ ean be seen as a funetion 

T : e ^ J (46) 

Considering the loop depleted at lines 5 to 9 in Algorithm]^ it ean be seen 
that every edge e e E{H) is eonverted in a pair of eomposite vertiees (u and 
v) and then represented as an edge on the adjaeeney matrix J(//). Therefore, 
if the eomposite vertiees u and v are adjaeent in MAG H, then a entry 1 is 
present at the interseetion of row Z)(u, t(//)) and eolumn of 
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J{H), indicating the corresponding adjacency in the matrix. Hence, the 
adjacency structure of the MAG H is preserved by the function T. 


Given the adjacency matrix J{H) and companion tuple we construct 
MAG H, which we then show to be isomorphic to the MAG H. We obtain 
A{H) from t{H) by constructing a list A(.^) with p = \ elements, in 

which every element i of this list is a set such that |A(.^)[/] | = t[z]. Without 
loss of generality, we can assume that the elements of each aspect A (.^)[z] 
are natural numbers ranging from 1 to t[z]. We then construct the edge set 
E{H), by starting with an empty set and then inserting an edge for each 
entry with value 1 in J{H). For constructing each of these edges, we take 
an entry of value 1, make r equal to its row number, and c equal to its 
column number. We then build the edge e = (Z)“^(r,T(//)),Z)“^(c,T(//)), 
which is a tuple of length 2p, where the first p entries correspond to the 
origin composite vertex and the last p entries correspond to the destination 
composite vertex of the edge. Note that function D~^ simply retrieves the 
original composite vertex entries from the row and column numbers of the 
adjacency matrix. Further, we construct the set V(.^) of composite vertices, 
which is the cartesian product of the sets in A(.^), so that 

p 

Y{H)= XMH)[nl (47) 

n=l 

where p = \t{H) \ is the number of aspects in the MAG H. 

We now show that the MAG H, constructed from J(//) and t(//), is isomor¬ 
phic to the original MAG H. Note that by construction of H we have that 

\E{H)\ = \E{H)\-, \A{H)\ = \A{H)\ = p; for I ^i^p,\A{H)[i\\ = |A(//)[z]|; 
|V(^)| = |V(//)|; and t(^) = x{H). 

Since |V(//)| = |V(.^)|, we know that there is a bijective function from 
Y{H) to Y{H). Further, we also have the bijective function D, which 
takes a composite vertex into a natural number, assigning a unique and dis¬ 
tinct natural number to each element of Y{H) and Y{H). Moreover, since 
x{H) = x{H) and by construction of D, we have that the range of D for 
Y{H) and V(.^) is the same, i.e. D{Y{H),x{H)) = D{Y{H),x{H)). From 
this, we conclude that, for every composite vertex u e Y{H), there is one 
unique composite vertex u eY{H) such that D(u, x{H)) = D(u, x{H)). We 
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( 48 ) 


thus define the bijeetive funetion 
/ : Y{H) Y{H) 

suehthatD(u, t(H)) =D(u, 

As the funetion / is bijeetive, for every edge e s H, we have an edge e = 
{f{Ko{e)), f{Kd,{e))) e E{H), and also, for every edge e e E{H), we have 
the eorresponding edge e = {f~^{no{e)),f~^{nd{e))) e E{H). This fulfils 
the eonditions for isomorphism between H and H. 

Sinee H is a quotient set under the MAG isomorphism relation and H is 
isomorphie to H, it follows that H and H eorrespond to the same element 
in H, making the funetion T bijeetive. Also, sinee eaeh entry with value 
1 in the adjaeeney matrix J{H) eorresponds to an edge in the MAG H, it 
follows that also preserves the MAGs adjaeeney strueture, establishing 
the isomorphism relation as desired. 

□ 

Theorem 2. Every function that can be obtained from a MAG to a given co¬ 
domain set can also be obtained from a matrix representation of the MAG. 

Proof Consider the diagram depieted in Figure In this figure, El is the set of 
all MAGs (up to isomorphism), J is the set of pairs of adjaeeney matriees and 
eompanion tuples (up to permutation), E is an arbitrary funetion from El to X, 
where X is a eodomain eonsistent with the definition of funetion E, and 1 is the 
identity funetion in X. Sinee the funetion E is arbitrary, it ean represent any 
funetion or algorithm, sueh as searehes or eentrality eomputations, whieh take 
MAGs to a result expeeted from this funetion. 


e 


T 


X 


X 


Figure 5: Commutative diagram. 


As both funetions T (Equation (46) in Theorem and 7 represent isomor¬ 
phisms, it follows that the depieted diagram eommutes, so that for every funetion 
F : El — > X there is a funetion F : J ^ X, whieh produees the same result. □ 
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As a consequence of Theorem it follows that, from the adjacency matrix 
and companion tuple of a MAG, one can obtain any possible outcome that can be 
obtained from a MAG or from any other representation equivalent to it, such as 
high order tensors, as those presented in recent related works p^[38l[^|. 


4.3 Degree 

The definition of degree in a traditional graph stems from the number of edges 
incident to a given vertex. This concept can be generalized for MAGs, so that 
degrees can be defined for composite vertices, sub-determinations, or elements 
of a given aspect. Further, since MAG edges are considered to be directed, the 
degrees are also divided into out-degree and in-degree. In this section, we present 
algorithms for calculating these distinct degree definitions. 

4.3.1 Degree of composite vertices 

The degree of composite vertices of a given MAG H can be obtained directly 
from its composite vertices representation, g{H). Since the composite vertices 
representation is a traditional directed graph isomorphic to the MAG H, it follows 
that the degree determination is done with the traditional algorithm for directed 
graphs with minor changes. For a given MAG H and its companion tuple x{H), 
the degrees of the composite vertices can be determined by Algorithm where 
D{Ko{e), t{H)) and D{nci{e), t(//)) stand for the numerical representation of the 
origin and destination composite vertices of edge e e E{H), as defined in Sec¬ 
tion |3^ 

Another way for calculating the degrees of the composite vertices is comput¬ 
ing it algebraically from the adjacency matrix of the MAG, as given by 


indegree = J(77)^ 1, 

(49) 

outdegree = J{H) 1. 

(50) 


Further, the total degree of the composite vertices can be obtained by summing up 
their indegrees and outdegrees. 

To determine the time complexity of Algorithm we consider that lines 2, 
4, and 5 have each time complexity G(|V(//) |), the determination of the compan¬ 
ion tuple at line 3 has complexity 0{p), where p is the number of aspects of the 
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input ■.H = {A,E) 
output: indegree, outdegrees 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


Degree(//) 

n-|V(H)| 

T ^ CompTuple{A[H)) //companion tuple of H, i.e. t{H) 

indegree •«— vector of n integers, all 0 
out degree •«— vector of n integers, all 0 
for each ee E{H) do 

o *— D[Ko{e),T) //numerical origin 

d <— D^TTd^e), T) //numerical destination 

indegree[d] <— indegree[d] + 1 
outdegree[o] <— outdegree[o] + 1 

end 

return indegree, outdegree 


Algorithm 8: Determination of the degree of composite vertices. 


MAG, so that p « |V(//)|. Finally, since the determination of the numerical rep¬ 
resentation of vertices has complexity 0{p), we have that the for loop initiated at 
line 6 has complexity 0{p *\E{El)\), so that the time complexity of Algorithm]^ 
is D(|V(//)| + p* \E{H)\). If we consider that in a given case the order of the 
MAG does not vary, so that p is a constant, then the algorithm’s time complexity 
isO(|V(H)| + |£(H)|). 

In the case of the example MAG T (Figure [^, whose companion tuple is 
T = (3,2,3), it can be seen that the composite vertex {2,Bus,tl) has outdegree 3 
and indegree 1, while the composite vertex {I,Subway,t2) has outdegree 2 and 
indegree 2. Since D((2,5M5,tl), t) = 2 and D{{1,Subway, t2),T) = 10, it follows 
that indegree[2] = 1, outdegree[2] = 3, indegree[lO] = 2 and outdegree[lO] = 2. 

4.3.2 Degree of sub-determined vertices 

We can determine the degree for sub-determined composite vertices in a similar 
way to the degree of composite vertices. Given a MAG H and a sub-determi¬ 
nation the degree of the sub-determined composite vertices can be obtained 
by Algorithm!^ where |V|^(//)| is the number of ^ sub-determined composite 
vertices on MAG H, S^ is the function that takes a composite vertex to its sub¬ 
determined form, and is the function that takes the sub-determined composite 


38 





vertex to its numerieal representation. It ean be seen that the time eomplexity of 
Algorithm]^ is the same as the time eomplexity of Algorithmic 


input : H = (A,E), and ^ 
output: indegree, outdegree 

1 SubDetDegree(//, 0 
|Vf(H)| 

CompTuple{A{H)) 

Tj ^ Ij(H) 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


n < 
T 


//companion tuple ofH, i.e. t{H) 
// sub-determined companion tuple 


indegree «— veetor of n integers, all 0 
outdegree «— veetor of n integers, all 0 
for each ee E{H) do 

o <— D[7ro{s),T^) //numerical sub-determined origin 

d D(7td{e),T^) //numerical sub-determined destination 

indegree[d] <— indegree[d] + 1 
outdegree[o] <— outdegree[o] + 1 

end 


13 return indegree, outdegree 

Algorithm 9: Sub-determined degree. 


It is important to note that two distinet eomposite vertiees may have the same 
sub-determined form. This happens when the two eomposite vertiees differ only 
on aspeets whieh are dropped by the sub-determination. In this ease, the degree of 
eaeh of these eomposite vertiees is summed for obtain the sub-determined degree. 
From this, it ean also be seen that some edges in the sub-determined form may 
beeome self-loops. The degrees ealeulated by Algorithm |C inelude the self-loop 
edges. This algorithm ean be modified to eount the self-loops separately, as shown 
in Algorithm [JC This algorithm is similar to Algorithm [C and has the same time 
eomplexity. 

The sub-determined eomposite vertiees degree ean also be determined alge- 
braieally with 


and 


indegree = 1, 

(51) 

outdegree = J{H) 1, 

(52) 
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input : H = {A,E), and ^ 

output: indegree, outdegree, self degree 

1 SubDetDegreeSepLoops(//, ^ 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


n <- 

r < 


|Vf(H)| 

CompTuple{A{H)) 

- hw 


//companion tuple ofH, i.e. t{H) 
// sub-determined companion tuple 


indegree •«— vector of n integers, all 0 
outdegree «— vector of n integers, all 0 
self degree «— vector of n integers, all 0 
for each e e E (H) do 


o ^ D{noie),Ti^) 
d ^D{Kd{e),T^) 

ltd ¥= o then 

indegree[d] <— 
outdegree\o] ^ 

end 

else 

I selfdegree[o] 

end 


// numerical sub-determined origin 
// numerical sub-determined destination 


indegree[d] + 1 
- outdegree[o] + 1 


selfdegree[o] + 1 


end 


19 return indegree, outdegree, self degree 


Algorithm 10: Sub-determined degree, separating self-loops 


where is the sub-determination matrix and 1 is the all Is column vector, 

both defined in Section 4.1 Note that the multiplication by adds the 

degrees of the composite vertices that are collapsed to the same sub-determined 
vertex. 

The degrees calculated by Equations ( |5T] ) and ( [5^ include the self-loop edges. 
To obtain the separate self-loop degrees, first note that 

i{H) = M^{H) i{H) 1„. (53) 


This follows from the fact that 


= (54) 

since M^(//)^ is a n x m rectangular matrix and has the property that each row 
has exactly one non-zero entry of value 1. 
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Furthermore, note that the matrix J{H) is the adjaeeney ma¬ 

trix of the sub-determined MAG Sinee the eomposite vertiees representation 
of a sub-determined MAG is a multigraph, eaeh non-zero entry shows the number 
of superposed edges in the sub-determination. Therefore, the main diagonal of 
J{H) has the self-loop degree of eaeh vertex. Henee, 


self degree = Diag{Mi^{H) i{H) (55) 


For example, eonsider the example MAG T (Figure]^ and the sub-determination 
= 01 12 defined in Seetion 4.1 We have that 


and 


indegree = M^,(r) J(T)^ 1 = 


outdegree = J(r) 1 = 


0 

7 

4 

4 

7 

0 

0 

7 

4 

4 

7 

0 


M{,(r)j(r) M{,(r4- 


0 0 0 0 0 0 

0 2 2 0 3 0 

0 2 2 0 0 0 

0 0 0 2 2 0 

0 3 0 2 2 0 

0 0 0 0 0 0 


self degree = Diag{M.i^f{T) J(r) M|^,(r)^) 


0 

2 

2 

2 

2 

0 


(56) 


(57) 


(58) 


(59) 


This means that, for instanee, the sub-determined eomposite vertex (2, Subway) 
has outdegree 7, indegree 7, and 2 self-loops. This sub-determination eorresponds 
to the aggregation of all 3 time instants, whieh means that the edges in whieh only 
the time instant ehanges beeome self-loops. These edges are shown in red (dotted) 
in Figure!^ Note that = (2,3), so that D{{2,Subway), (2,3)) = 2, making it 
eorrespona to the seeond element of the degree eolumn veetor. 
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4.3.3 Single aspect degree 


The single aspect degree is a particular case of sub-determined degree in which 
the sub-determination applied is such that only a single aspect remains. Therefore, 
the determination of single aspect degrees is done in the same way presented in 
Section l4.3.2[ 

We, however, present an additional example illustrating the time instant de¬ 
gree, which is obtained by the sub-determination = IOO 2 defined in Sec¬ 
tion |4.1[ This sub-determination has only the third aspect of the MAG T (Fig¬ 
ure [^, which corresponds to the three time instants present on MAG T. In this 
case, we have that 


indegree = M|^ 7 ’(r) J(T)^ 1 = 


outdegree = M^ 7 '(r) J(r) 1 = 


Mfr(r) W - 


2 
10 
10 

10 
10 
2 

8 0 
2 8 
0 0 2 


2 

0 


and 


self degree = Diag{fAi^j{T) J(r) = 


2 

2 

2 


(60) 


(61) 


(62) 


(63) 


Therefore, we have that = (3), so that T>((ti), (3)) = 1, £>((^ 2 ), (3)) = 2, and 
Z)((t 3 ), (3)) = 3. Considering the composite vertices representation of MAG T, 
depicted in Figure it can be seen that each time instant has 2 self-loop edges (in 


blue-dashed), which is consistent with Equation (63). Further, there are 8 edges 
from ti to t 2 (in red-dotted and black) and 8 edges from t 2 to . This is consistent 


with the adjacency matrix shown in Equation (62). Further, the indegrees and 


outdegrees of each time instant are consistent with Equations ([60]) and ([M]). 


4.4 Breadth-First Search (BFS) 

The Breadth-Eirst Search (BES) is an important graph algorithm that can be seen 
as a primitive for building many other algorithms [[^ . The goal of this section is 
to illustrate how the BES algorithm can be adapted for being used in MAGs, both 
in its full composite vertices representation and in its sub-determined forms. In 
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the not sub-determined form, the adaptation is very simple, sinee the eomposite 
vertiees representation of a MAG is a direeted graph. In this ease, all that is needed 
is to eonvert the eomposite vertiees representation from its tuple to numerieal 
form, and then apply the traditional BPS algorithm. The adaptation to the sub¬ 
determined forms also does not require major ehanges on the algorithm. As with 
many graph algorithms, BPS ean be expressed in eombinational or in algebraie 
forms, whieh are presented in the following related subseetions. 


4.4.1 BFS for composite vertices 


The non sub-determined BPS in its eombinational form is eonstrueted direetly 
upon the MAG’s adjaeeney matrix, J{H). 

Considering Algorithm and the standard form of the BPS algorithm en- 
eountered in [[^, it ean be seen that the differenee is that the starting eomposite 
vertex s has to be transformed from its tuple representation to its numerieal repre¬ 
sentation, as shown in lines 8,9 and 10 of Algorithm[^ Therefore, from the anal¬ 
ysis provided in Q, we ean eonelude that the time eomplexity of Algorithm [TT] is 


0(|V(H)| + |ii:(H)|). 

BPS is also elosely related to matrix multiplieation. This stems from the well- 
known property of the powers of the adjaeeney matrix, in whieh the (/, j) entry 
of the n-th power of the adjaeeney matrix shows the number of existing walks of 
length n from vertex i to vertex j [jM|. Prom this, we eould think that for a given 
MAG H, the series 


00 

B = YjJiHY = I + J{H) + J{Hf + J(//)3 -f J(//)4 + ... (64) 

i=0 


would produee a matrix B, sueh that the entry B,- ^ indieates the number of walks 
of any length from vertex i to vertex j. This is indeed the ease when H happens to 
be an aeyelie MAG, making J{H) a nilpotent matrix. 

The existenee of eyeles in H makes that, for some vertiees, there will exist 
walks of arbitrary length eonneeting them (namely, the eyeles), making the series 
of Equation (|6^ divergent. However, sinee the objeetive is not to know the num¬ 
ber of walks between eaeh pair of vertiees, but simply to know whieh vertiees are 
reaehable from eaeh other (i.e. there is at least a path between them), this teehni- 
eal problem ean be solved by multiplying the adjaeeney matrix J{H) by a sealar 
pH, sueh that 
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input : andse V(//) 

output: vertices, distance, pred 

1 BFS(J(//),t(//),s) 

2 n^|V(//)| 

3 vertices *— vector of n integers, all 0 

4 distance <— vector of n integers, all oo 

5 pred <— vector of n integers, all Nil 

6 color <— vector of n integers, all 0 

7 2 «— empty queue 

8 vertices\D{s, Tau) — 1] <— 1 

9 distance\D(s,x(H))\0 

10 Enqueue{Q,D{s,x{H))) 

11 while Q not empty do 

12 M«— head\Q] 

13 for each v successor of u do 

14 if co/or[v] = 0 then 

15 I co/or [v] 1 

16 vertices\v^ <— 1 

17 distance\v^ <— distance\u\ + 1 

18 pred\v^ <— u 

19 Enqueue{Q,v) 

20 end 

21 end 

22 Dequeue (Q) 

23 color[u] <— 2 

24 end 

25 return vertices, distance, pred 

Algorithm 11: BFS for composite vertices. 


where p(J(//)) is the spectral radius of the matrix This leads to the matrix 

Jp(//)=PH J(//), (66) 


so that the spectral radius of the matrix ip{H) < 1. This results that Equation (64) 
constructed with the matrix ip{El) converges. Since the convergence of the series 
is assured. Equation ( [64] ) can be re-expressed as 

B = (/-Jp(//))-i. (67) 
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The matrix B defined in Equation ( [FT] ) has the property that, for any given 
eomposite vertex v e Y{H), the row D{y) of B has non-zero entries in every eol- 
umn that eorresponds to a eomposite vertex u e Y{H), sueh that u is reaehable 
from v. Henee, for a given eomposite vertex v, the row D{y) eorresponds to the 
result of a BPS started at that eomposite vertex. Although the matrix B earries the 
BPS of all eomposite vertiees of the MAG H, it is important to note that this ma¬ 
trix may not be sparse, whieh for large MAGs ean lead to diffieulties in memory 
alloeation. In order to avoid sueh diffieulties, it is also possible to express a BPS 
for a single eomposite vertex v as 

B = Tyl + ryJp[H) + ryJp[H)^ + f\3p+ ryJp(H)‘^ -I- ..., (68) 


where is the row veetor with n entries for whieh all entries exeept D(v, t) are 0 
and the entry D(v, t) is 1. 

Considering the example MAG T, shown in Pigure the result of the BPS 
using Algorithm 11 for the eomposite vertex {2,BusJi), whose numerieal repre¬ 
sentation is (3,2,3)) = 2, is 


vertices = [2,5,8,9,10,11,14,15,16,17] (69) 

distances = [oo,0,oo,oo, l,oo,oo, 1,1,2,2,oo,oo,2,2,3,3,oo] 

pred = [MZ,MZ,M/,M/,2,M/,M/,2,2,5,5,MZ,M/,8,8,10,10,AZZ], 


where the list vertices shows the eomposite vertiees aeeessible from {2,Busd\), 
whieh in this example represent all loeations, transit modals, and time instants 
reaehable from this initial point. The list distances earries the distanees in hops 
from the initial eomposite vertex (2,Busd\) to all possible destinations (with oo 
meaning that a destination is not reaehable). The list pred shows the predeeessors 
of eaeh eomposite vertex, making possible to eonstruet a BPS tree. 


4.4.2 Sub-determined BFS 

It is possible to obtain a sub-determined form of the BPS algorithm for MAGs. 
It is important, however, to realize that this sub-determined BPS algorithm is not 
equivalent to applying the BPS algorithms presented in Seetion 4.4.1 to a sub¬ 
determined MAG. A sub-determination is a generalization of the idea of aggre¬ 
gating multilayer and time-varying graphs, as shown in Seetion |2.2[ As with the 
aggregation proeess, the sub-determination of a MAG ean eause the presenee of 
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(a) MAG R (b) Sub-determined form 

Figure 6: MAG R and its sub-determined form. 


paths and walks on the sub-determined MAG that do not aetually exist on the orig¬ 
inal MAG. To illustrate this, we present Figures and whieh show a small 
two aspeets MAG and its sub-determined form, obtained by the sub-determination 
^R = 0 12 . First, note that, in the MAG R shown in Figurethere is no path orig¬ 
inating from the eomposite vertiees (1,1) or (1,2) to the eomposite vertiees (3,1) 
or (3,2). Nevertheless, in Figure]^ there is a path eonneeting the sub-determined 
vertex (1) to the sub-determined vertex (3), even though sueh eonneetion is not 
possible on the original MAG shown in Figure Therefore, in order to obtain 
the proper result, the sub-determined BFS should not be evaluated direetly using 
the sub-determined MAG. 

Sueh a ease ean be seen algebraieally by noting that given a MAG H and a 
sub-determination in general 


00 


00 


M?(H) X;Jp(H)M M{(Hy/2(M5(H)Jp(H)Mf(ff)'')'. (70) 


Ki=0 


i=0 


To see that the Inequality (70) holds, note that an arbitrary power of the matrix 
M^{H) Jp(//) M^(//)^ is givenby 

(M^(//)Jp(//)M^(//)^)" = M^(//)Jp(//)M^(//)^M^(//)Jp(//)M^(//)^.., 


V 

n times 


where (M^(//) Jp{H) is multiplied n times. Note, however, that 

M^{Hf M^{H) 


(71) 


(72) 
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since M|r(//) e is a retangular matrix and m < n, so that the rank of the 

matrix is less or equal to m, while the rank of the identity I„ is 

n> m. Sinee Inequality ( [72] ) holds, so does the Inequality ( |70| ). 

Here, the left hand side of the Inequality fTO] ) eorresponds to the sub-determina¬ 
tion of the BPS ealeulated for the MAG H, while the right hand side eorresponds 
to the BPS ealeulated for the sub-determined MAG 

In the ease of the MAG R, shown in Pigurej^ we have that the sub-determina¬ 
tion is given by = 0 12 and the adjaeeney and sub-determination matriees are 


JW = 


0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


(73) 




10 0 10 0 
0 10 0 10 
0 0 1 0 0 1 


and 


= 


0 10 
0 0 1 
0 0 0 


Therefore, we have that 


00 




yi=0 


while 


and 


00 




i=0 


3 2 0 
0 3 2 
0 0 3 


2 2 2 
0 2 2 
0 0 2 




1 1 1 
0 1 1 
0 0 1 


(74) 


(75) 


(76) 


(77) 


(78) 


Remembering that the entries of the matriees in Equations ( |7^ , ( [77] ), and ( [78] ) are 
to be eonsidered only as zero or non-zero, it ean be seen that the matrix at Equa¬ 
tion (76) has a 0 at entry (1,3), while the matriees at Equations ( [77j ) and (78) have 
a non-zero entry at this same position. This illustrates the situation in whieh a BES 
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is done on the sub-determined (aggregated) MAG, as in Equations ( |77| ) and ( [78] ), 
i.e. paths that are not present on the original MAG ean appear on the sub-determi¬ 
ned form, potentially altering the results obtained by algorithms applied to it. 

For instance, considering the MAG T, depicted in Figure]^ for a sub-determi¬ 
nation = 01 12 , which drops the time aspect, and considering pn = 0.5 so that 
Jp(r) = 0.5 J(r), we have that 


Mc,(r) 



Mf,(?■)'■ 


3 

0 

0 

0 

0 

0 

0 

7.8 

2.2 

1.3 

5.2 

0 

0 

2.2 

4.6 

0.2 

1.3 

0 

0 

1.3 

0.2 

4.6 

2.2 

0 

0 

5.2 

1.3 

2.2 

7.8 

0 

0 

0 

0 

0 

0 

3 


(79) 


Algorithm[^shows a combinational version of the sub-determined BFS. This 
procedure ensures that only paths present on the original MAG are considered on 
the sub-determined BFS. The sub-determination of the results obtained from the 
BFS is done in the internal if, comprising lines 21 to 25 of Algorithm [T^ 

After applying Algorithm to the MAG R with initial vertex ^ = 1 and sub- 
detemination = 0 12 , the obtained result is 


vertices = [1,2] (80) 

distances = [0,l,oo] 
pred = [A//, I,Nil], 


which is consistent with the result obtained by Equation ( fT^ . 

Further, applying Algorithm to MAG T, shown in Figure with start¬ 
ing composite vertex 5 = [2,Bus) and applying the sub-determination = 01 12 , 
which drops the time aspect, the obtained result is 

vertices = [2,5,3,4] (81) 

distances = [ 00 , 0 ,1 ,2, 1 , 00 ] 

pred = {Nil,Nil,2,5,2,Nil]. 


Considering that Xq = (3,2) and 2 = D{{2,Bus), (3,2)), 5 = D{{2, Subway), (3,2)), 
3 = D{{3,Bus), (3,2)), and 1 = D{{1,Subway), (3,2)), this means that disregard¬ 
ing time, starting from [2,Bus) it is possible to reach {2, Subway) in 1 step. 
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input : J{H), t{H), ^ and s e V^(//) 
output: vertices, distance, pred 
1 BFS-Sub(J(//),T(//),C,s) 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


\Y{H)\ nS- 




// ^ sub-determined companion tuple 


vertices <— vector of nS integers, all 0 
distance <— vector of nS integers, all oo 
pred <— vector of nS integers, all Nil 
colors <— vector of nS integers, all 0 
color <— vector of n integers, all 0 
Q <— empty queue 

for every \eY{H) where = D{s,z^{H)) do 

Enqueue{Q,D{\, t(//))) 
color[D{\,T{H))] ^ 1 

end 

vertices[D{s,Tz) — 1] <— 1 
distance[D{s,Ti^)] <— 0 
while Q not empty do 
u <— head[Q] 

for each v successor of u do 
if color\v\ = 0 then 
color\v^ <— 1 
Enqueue{Q,v) 
if colorS[D{v, T^)] =0 then 
colorS[D{v,T^)] <— 1 
vertices[D{v,T^)] <— 1 
distance[D{v, T^)] <— distance[D{u, T^)] 
pred[D{v, T^)] <— D{u, T^) 
end 
end 
end 

color[u] <— 2 
Dequeue{Q) 

end 


33 return vertices, distance, pred 


Algorithm 12: Sub-determined BPS. 
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{I,Subway) in 2 steps, and {3,Bus) in 1 step. It is not possible to reaeh (l^Bus) 
beeause there is no bus stop at loeation 1, neither (3,Subway) beeause there is 
no subway station at loeation 3. From the predeeessor list (pred) it is possible 
to build a BFS tree, where {2,Bus) is the root, {2,Subway) and {3,Bus) are ehil- 
dren of {2,Bus), and [I,Subway) is a ehild of (2,Subway). Note that {I,Subway) 
and {3,Bus) are leaves. It ean be seen that the result obtained in Equation (79) is 
eonsistent with the results obtained by Algorithm [T^ Comparing Algorithm [T^ to 
Algorithmic it ean be seen that the main differenee is the additional for loop at 
line 12 of Algorithm [C Sinee the time eomplexity of this loop is 0(|V(//)|), we 
then eonelude that the time eomplexity of Algorithmic is 0(|V(//)| + \E{H)\). 


4.4.3 Single aspect BFS 

The single aspeet BFS is a speeial ease of the sub-determined BFS. As sueh, it 
is evaluated using the same algorithms presented in Seetion |4.4.2| for the sub¬ 
determined ease. 

Considering the example MAG T (Figure |^, a sub-determination = 0 0 12 , 
whieh drops the time and transit mode aspeets (thus leaving only the loeations 
aspeet), and making pH = 0.5 so that ip{T) = 0.5 J(r), we have that 




w=0 


7.6 

3.5 

0.2 ■ 

3.5 

26 

3.5 

0.2 

3.5 

7.6 


(82) 


indieating that disregarding the aspeets of time instants and transit modes, all 
loeations ean be reaehed from any loeation. 

Applying Algorithm |C to the MAG T, with starting eomposite vertex 5 = (1) 
and employing the sub-determination t^i= OOI 2 , whieh drops the aspeets of the 
transit mode and time instants, the obtained result is 


vertices = [1,2,3] (83) 

distances = [0,1,2] 
pred = {Nil, 1,2], 


whieh is eonsistent with the result obtained by Equation (82). 
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4.5 Depth-First Search (DFS) 

In this section, we show the adaptation of the Depth-First Search (DFS) algo¬ 
rithm for use with MAGs. The DFS algorithm exposes many properties of the 
MAG structure and can be used as a primitive for the construction of many other 
algorithms Q. We present DFS algorithms for both the full composite vertices 
representation of the MAG as well as for the sub-determined form. We remark 
that in the sub-determined algorithm the full information of the MAG is used, in 
the sense of preventing the use of paths that may exist in the sub-determined form 
of the MAG, while not actually existing in the original MAG. 

4.5.1 DFS for composite vertices 

The composite vertices implementation is constructed using the MAG’s adjacency 
matrix J(//) and companion tuple t(//). The implementation shown is very simi¬ 
lar to the traditional implementation presented in Q, which is expected since the 
composite vertices representation of the MAG is indeed a directed graph, so that 
the original algorithm applies. 

The proposed implementation can be seen in Algorithm is similar to the 
original implementation. Therefore, considering the analysis provided in Q, we 
conclude that the time complexity of Algorithm ^is D( | V(//) | -f jF (H) \). 

When applied to MAG T, shown in Figure [3 the DFS algorithm generates the 
result 


d = [0,2,22,24,3,26,28,13,19,4,12,30,32,8,14,5,7,34] (84) 

/= [1,21,23,25,18,27,29,16,20,11,17,31,33,9,15,6,10,35] 
pred = [Nil,Nil,Nil,Nil, 2,Nil,Nil, 11,2,5,5,Nil,Nil, 17,8,10,10, W7], 


where the list d carries the discovery time of each composite vertex, the list / the 
respective finish time of each composite vertex, and pred the predecessor list of 
each composite vertex. 

4.5.2 Sub-determined DFS 

The sub-determined DFS algorithm is presented in Algorithm [T^ and is similar to 
the non sub-determined one. The main differences are at the Procedure Visit-DFS- 
Sub and the call to a sub-determined BFS at line 15 of the DFS-Sub function. This 
version for a sub-determined BFS is considered in order to determine reachability 
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input :J{H),x{H) 

output: discTime, finTime, pred 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


DFS(J(//), t{H)) 
n^\V{H)\ 
for u = I to n do 


color\u\ 0 

// set all vertices to white 

discTime\u\ * -1 

// set discovery times to nil 

finTime[u] < -1 

//setfinish times to nil 

pred[u] * -1 

// set predecessors to nil 

end 


time <— 0 


for u = \ to n do 



if color[u] = 0 then 
I DFS-Visit{u) 

end 

end 

return discTime, finTime, pred 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


Procedure DFS-Visit(u) 
color[u\ 1 
discTime[u\ *—time 
time <— time + 1 
for each v successor ofu do 
if color{v\ = 0 then 
pred[v\ <— u 
DFS-Visit{v) 
end 
end 

color[u] <— 2 
finTime\u\ <— time 
time <— + 1 


// set vertex u to gray 


// set vertex u to black 


Algorithm 13: DFS for composite vertices. 


of sub-determined vertices from the root of each sub-determined DFS tree. This 
is necessary to prevent including vertices not reachable from the tree root in the 
non sub-determined MAG into the DFS trees constructed by Procedure Visit-DFS- 
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Sub. An example of this is provided in Equation ( [8^ . The differenee in Proeedure 
Visit-DFS-Sub is that in addition to the root vertex for the DFS tree it also reeeives 
the reaehability veetor produeed by the BFS. This reaehability veetor has one entry 
for eaeh sub-determined vertex. This entry has value 1 when eorresponding to a 
reaehable vertex, while entries eorresponding to unreaehable vertiees earry value 
0 . 

In order to determine the time eomplexity of Algorithm [T^ we eonsider that 
the sub-determined BFS exeeuted at line 15 of Funetion DFS-Sub is done onee 
for the root vertex of eaeh sub-determined DFS tree. Sinee it is exeeuted only 
onee for eaeh DFS tree, we eonelude that the total time expended in the sub¬ 
determined BFS algorithm is D(|V(//)| -I- \E{H)\). Sinee the reaehability eheek 
ineluded in Funetion Visit-DFS-Sub is done by verifying the eontent of one entry 
of the reaehability veetor, it is done in 0{l) and therefore does not affeet the 
overall time eomplexity of the Visit-DFS-Sub Funetion. Therefore, sinee the DFS 
is run upon the sub-determined MAG, it follows that the time eomplexity of doing 
the DFS part of the Algorithm is D(|V|j(//)| -t-1£’^(//)|). Sinee |V^(//)| < |V(//)| 
and \Ei^{H)\ < \E{H)\, we eonelude that the time eomplexity is dominated by the 
BFS used for the reaehability determination, making the overall time eomplexity 
of Algorithm [14] to be D(|V(//)| -I- \E{H)\). 

When applying the sub-determined DFS algorithm to the example MAG T 
shown in Figure [^ with a sub-determination Q = OII 2 , whieh drops the time 
aspeet, the obtained result is 

[0,2,3,6,5,10] (85) 

/= [1,9,4,7,8,11] 
pred = [Nil,Nil,2,5,2,Nil], 


where the list d earries the diseovery time of eaeh sub-determined eomposite ver¬ 
tex, the list / its finish time and pred its predeeessor. 

Considering the MAG R shown in Figure 6a with a sub-determination Q; = 
0 12 , the result obtained by Algorithm [T4| is 


[0,1,4] (86) 

/= [3,2,5] 
pred = [Nil, 1,V//]. 


It ean be seen that even though in the MAG R sub-determined by = 0 12 (see 
Figure |6^ there is a path from vertex 1 to 3, vertex 3 is not in the same DFS 
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inpul 




output: discTime, finTime, pred 


1 

E 

FS-Sub(J(//), t{H),0 


2 


T 


// ^ sub-determined companion tuple 

3 


<— SubDetMatrix{H, Q 


4 


J 

f = Mf i{H) M[ 

// Sub-determined adjacency matrix 

5 


n 


//number of sub-determined vertices 

6 


for u = 1 to n do 


7 



color\u) <— 0 

// set all vertices to white 

8 



discTime\u) <-1 

// set discovery times to nil 

9 



finTime\u) <-1 

//setfinish times to nil 

10 



pred[u) <-1 

// set predecessors to nil 

11 


end 


12 


time <— 0 


13 


for u = 1 to n do 


14 



if color[u] = 0 then 


15 



vertices = BFS-Sub{J{H), 


16 



DFS-Visit-Sub{u,vertices) 


17 



end 


18 


end 


19 

rtinrn discTime, finTime, pred 


1 

Procedure DFS-Visit-Sub(u, vertices) 


2 


color[u] <— 1 

// set vertex u to gray 

3 


discTime[u\ <— time 


4 


time <— time + 1 


5 


for each v successor of u do 


6 



if color[v] = 0 and vertices[v] ^ 0 then 

7 



pred[v] <— u 


8 



DFS-Visit-Sub{v, vertices) 


9 



end 


10 


end 


11 


color\u) <— 2 

// set vertex u to black 

12 


finTime\u\ <— time 


13 


time <— time + 1 



Algorithm 14: Sub-determined DFS. 


tree as vertiees 1 and 2, even with the DFS starting at vertex 1, as ean be seen in 
J[0]. This oeeurs beeause in MAG R (with no sub-determination) there is no path 
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connecting the composite vertex 1 to the composite vertex 3. 


4.6 Single aspect DFS 

The single aspect DFS is a special case of the sub-determined DFS. As such, it 
is evaluated using the same algorithms presented for the sub-determined case in 
Section l4.5.21 

Applying Algorithm[^to MAG T (Figure]^, and employing the sub-determi¬ 
nation = 0 0 12 , which drops the aspects of transit modes and time instants, the 
obtained result is 


[0,1,2] (87) 

/= [5,4,3] 
pred = [A//, 1,2], 


where the list d carries the discovery time of each sub-determined composite ver¬ 
tex, the list / its finish time, and pred its predecessor. 


5 Final remarks 


In this paper, we have presented the algebraic representation and basic algorithms 
of MultiAspect Graphs (MAGs). The key contribution has been to show that mod¬ 


els based on the MAG abstraction (formally defined in [ 19|) can be represented 
by a matrix and a companion tuple. Furthermore, we have also shown that any 
possible MAG function (algorithm) can be obtained from this matrix-based rep¬ 
resentation. This is an important theoretical result because it paves the way for 
adapting well-known graph algorithms for application in MAGs. In this sense, 
we have presented the adaptation for the MAG context of basic graph algorithms, 
such as computing degree, BFS, and DFS. In particular, we have also presented the 
sub-determined versions of the same basic algorithms, showing that such versions 
disregard spurious paths that usually result from the sub-determination process, 
thus avoiding the pollution of the results with the consideration of such paths. 

As future work, we intend to build upon the results here obtained for the alge¬ 
braic representation and basic algorithms of MAGs to analyze MAG properties, 
such as the centrality of edges, composite vertices, and aspects. We also intend to 
consider the dynamics encountered in these properties in the cases where one of 
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the MAG aspects represents time. Finally, we are also targeting the application of 
the MAG concept for the better understanding, modeling, and analysis of different 
complex networked systems found in real-world applications. 
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